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Se presenta la validación de una metodología de agre- 
gación, propuesta para obtener conductividades hidráu- 
licas equivalentes de un medio homogéneo en suelos 
heterogéneos no saturados. Para ello se propone, en una 
primera etapa, una modelación numérica del flujo en la zona 
no saturada, en el dominio heterogéneo discretizado con 
alto grado de detalle y parametrizado mediante funciones 
de pedotransferencia adecuadas. En una segunda etapa se 
calibran los parámetros equivalentes del modelo de frente 
húmedo con los resultados alcanzados en la primera etapa, 
comparándolos con los propuestos en la metodología de 
agregación. El procedimiento se aplica a 68 perfiles típicos de 
suelos que caracterizan el sector meridional de la provincia 
de Santa Fe (Argentina). Los resultados alcanzados validan 
la metodología de agregación propuesta para un numeroso 
texturales, 
hidráulicas y contenidos de humedad antecedente. Para 


muestreo de composiciones propiedades 


las series de suelos analizadas, serían aceptables las 
estimaciones de conductividades hidráulicas saturadas 
mediante ponderación hasta saturaciones efectivas del 
orden del 50%, con un rango de errores acotados al 10%. 


Palabras clave: suelos heterogéneos no saturados, conduc- 
tividad hidráulica, metodología de agregación. 


Introducción 


Para predecir el movimiento de agua y solutos 
en la zona no saturada (ZNS) se emplean 
habitualmente modelos numéricos. Los 
resultados de sus aplicaciones sólo pueden 
ser confiables si las propiedades del medio 
poroso son conocidas. Esto representa la 


principal limitación a la tecnología actual de 
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This study presents the validation of an aggregation methodology 
proposed to obtain hydraulic conductivity equivalents for a 
homogeneous porous media in unsaturated heterogeneous soils. 
The first step in this proposal is to numerically model flow in the 
unsaturated zone of a heterogeneous area, with a highly detailed 
discretization and parameterization based on suitable pedotransfer 
functions. The second stage involves calibrating the equivalent 
parameters from the wet front model based on the results from the 
first phase and comparing with those proposed by the aggregation 
methodology. The procedure was applied to 68 soil profiles 
characteristic of the southern province of Santa Fe (Argentina). 
The results obtained validate the aggregation methodology proposed 
for a large sampling of textural compositions, hydraulic properties 
and antecedent moisture contents. For the series of soils analyzed, 
estimates of saturated hydraulic conductivity would be acceptable 
with weighting up to effective saturations of 50%, with a margin of 
error limited to 10%. 
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modelos complejos que se encuentra en un 
estado avanzado con respecto al conocimiento 
de los parámetros que engloban dichos 
modelos. En trabajos de calibración de modelos 
hidrológicos continuos se ha demostrado la 
alta sensibilidad que presentan los procesos 
subterráneos —observables en el sistema 
mediante los niveles freáticos— y los procesos 
—observables a 


hidrológicos superficiales 
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través de volúmenes de escurrimiento— frente 
a los parámetros hidráulicos que caracterizan 
la zona no saturada (Zimmermann, 2003). 
También estos parámetros inciden en la 
estimación de procesos de erosión y deposición 
de sedimentos a escala de cuenca (Basile, 
Riccardi, Zimmermann, éz Stenta, 2010). 

Para caracterizar el medio no saturado, 
se requiere del conocimiento de las curvas 
de conductividad K(8), y retención wyp(0) o 
bien, difusividad D(0), siendo O el contenido 
volumétrico de humedad, K(0) la conducti- 
vidad hidráulica no saturada, y(8) el potencial 
mátrico del suelo y D(0) el coeficiente de 
difusividad. Se ha diseñado un gran número 
de ensayos de laboratorio y métodos de 
campo para medir las funciones hidráulicas 
del suelo, pero la mayoría de los métodos son 
relativamente costosos y difíciles de llevar a 
cabo. 

Existen métodos indirectos para estimar las 
propiedades hidráulicas de suelos, que varían 
en términos de la metodología y complejidad, 
pudiendo distinguirse tres grupos principales: 
los métodos basados en la distribución del 
tamaño de poros, los métodos inversos y las 
funciones de pedotransferencia. 

Las funciones de pedotranferencia (EPT) 
ofrecen una metodología para aproximar 
las propiedades hidráulicas, empleando la 
composición textural del suelo, además de otra 
información disponible (p. ej., la distribución 
del tamaño de partículas, densidad aparente, 
contenido de materia orgánica, etc.). Existen 
numerosas funciones de pedotransferencia, 
explicitadas como ecuaciones y/o imple- 
mentadas como códigos informáticos, que 
han sido ajustadas a diferentes bases de 
información edafológica. En trabajos previos 
se han utilizado FPT para estimar propiedades 
hidráulicas de suelos limosos representativos 
de un sector de la llanura argentina, partiendo 
deinformación granulométrica y de parámetros 
físico-químicos disponibles en las cartas de 
suelos elaborados por el Instituto Nacional de 
Tecnología Agropecuaria (INTA). Se emplearon 
FPT incluidas en los códigos SoilPar y Rosetta, 


junto con ecuaciones de regresión entre 
parámetros hidráulicos y parámetros físicos 
(Zimmermann éz Basile, 2007). 

A partir de comparaciones entre parámetros 
obtenidos de FPT 
y los ajustados en modelos matemáticos 


mediante aplicaciones 


hidrológicos, pudo concluirse que de los 
códigos empleados, el programa Rosetta 
(Schaap, Leij, £ Van Genuchten, 1999), fue el 
que mejor se aproximó para la zona de estudio 
(Zimmermann éz Basile, 2011). 

La siguiente fase del trabajo consistió 
en la caracterización hidráulica de suelos 
de composición heterogénea. Los suelos no 
constituyen medios homogéneos, sino que 
los perfiles se presentan estratificados en 
horizontes y además se combinan formando 
asociaciones y consociaciones. Por ende, con 
el propósito de simplificar el abordaje del 
modelado, fue propuesta una metodología 
la estimación de 
parámetros representativos, desde un punto 


de ponderación para 


de vista hidráulico, equivalentes a un medio 
homogéneo (Zimmermann et al., 2008). Esto 
se hizo con el propósito de utilizar dichos 
parámetros en modelos simplificados de 
infiltración, tales como el propuesto por Green 
y Ampt (frente húmedo). 

En este trabajo se pretende validar esta 
metodología de ponderación, contrastando 
los resultados obtenidos con un modelo 
matemático basado en la resolución de la 
ecuación de Richards, en donde se contempla 
la heterogeneidad de cada perfil de suelo, el 
cual se describe a continuación. 


Modelo numérico empleado para la 
simulación del flujo en la ZNS 


El modelo de flujos empleado para las 
denominado SOLUM, está 
basado en la ecuación de Richards y se desa- 


simulaciones, 


rrolla específicamente para este trabajo, con el 
fin de estimar la redistribución de humedad 
en la zona no saturada (ZNS) y los volúmenes 
de agua intercambiados con la atmósfera y el 
acuífero. El mismo forma parte de un modelo 
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más general, desarrollado para simular el 
movimiento de flujos superficiales en la zona 
saturada y la no saturada, en áreas de llanura, 
denominado SHALL (Simulación Hidrológica 
de Áreas de Llanura) (Zimmermann é Riccardi, 
2000). En el modelo, la resolución de la ecuación 
de flujo en la ZNS se realiza en términos del 
contenido volumétrico de humedad 0, para 
conocer una de las variables del balance 
hidrológico de manera directa. Solamente está 
contemplado el flujo en la dirección vertical. 
La ecuación de movimiento, en términos del 
contenido volumétrico de humedad y con una 
analogía a la ley de difusión de Fick, es: 


(1) 


donde q es la velocidad de Darcy del flujo 
no saturado; 6, el contenido volumétrico de 
humedad; z, la coordenada vertical; zr, la 
profundidad radicular; K(6), la conductividad 
hidráulica no saturada; D(0), el coeficiente 
de difusividad del medio poroso, producto 
de la conductividad no saturada K(8) por el 
gradiente de la curva de retención (dW /d0), 
y S(8) un término sumidero que, en el caso 
de suelos cultivados, representa la tasa de 
extracción de agua por las raíces de las plantas. 
El modelo prevé incluir la heterogeneidad del 
suelo, otorgando distintos parámetros que 
regulan las curvas de retención y conductividad 
hidráulica para los diferentes estratos del suelo. 
La aplicación que se le dio al modelo en este 
trabajo consistió en la simulación del proceso 
de infiltración en un evento lluvioso, por lo que 
el término sumidero no fue considerado. En la 
ecuación, un valor negativo de q denota flujo 
hacia abajo (dado que q es positivo cuando el 
flujo es en la dirección z positiva hacia arriba). 
La ecuación de continuidad establece para el 
flujo vertical que: 


%__%9 
ot 0Z 


2) 


El medio poroso no saturado se representa 
en forma discreta mediante un conjunto de 


celdas que se extienden verticalmente desde 
la superficie hasta el nivel freático. Para la 
resolución de las ecuaciones (1) y (2) se ha 
propuesto un esquema numérico de tipo 
explícito en diferencias finitas. El esquema 
es centrado en el espacio y progresivo en el 
tiempo. En los bordes de celdas se evalúan 
los flujos de intercambio y en los centros de 
celdas se estima el contenido de humedad. Con 
este esquema de discretización (figura 1), la 
ecuación de momentum resulta: 


dol EZ ao) 


0/+1/2 + 0-1/2 
HE 65) 


n_ 

0j = 

La ecuación de continuidad resuelta para el 
término en el nivel n+1 es: 
n+1 At 


n n 
Oj+1/2 = Ojr1/2 — laj=8ja) 


e (4) 


donde j y n representan los índices de 
discretización espacial y temporal, respecti- 
vamente. La secuencia de cálculo es la 
siguiente: a) partir de una condición inicial 
de contenidos de humedad Or, en el perfil; b) 


resolver la ecuación de momentum (3) para el 


Az 
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Figura 1. Discretización de la ZNS en celdas verticales. 
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espesor de la ZNS, determinando los flujos 
de intercambio entre celdas; c) resolver la 
ecuación de continuidad (4), determinando 
los contenidos de humedad del perfil en el 
siguiente paso de tiempo, y d) retomar el paso 
(b) hasta alcanzar el tiempo de finalización de la 
simulación. El código de resolución fue escrito 
enteramente en lenguaje FORTRAN. Para el 
esquema propuesto, se ha realizado un análisis 
de estabilidad y convergencia, contrastando sus 
resultados con problemas de solución analítica 
conocida y con otras técnicas y esquemas 
numéricos, validando su consistencia en 
todas las situaciones (Zimmermann, 1998). 
De manera complementaria, se ha realizado 
un análisis de estabilidad, observando el 
comportamiento de números adimensionales 
(Courant € Péclet) en situaciones críticas de 
humedad, programando alertas en el código 
fuente para situaciones donde los números 
escapan del entorno de estabilidad numérica. 
Como condiciones de contorno, se plantea 
la resolución de las ecuaciones anteriores, 
con algunas de las incógnitas conocidas. La 
imbibición como condición de borde superior 
fue resuelta utilizando el esquema numérico 
como un modelo de infiltración. En la celda 
que constituye el contorno superior se debe 
cumplir con la siguiente condición: 


n 


ama o foja) fpa);-0" 6 


donde ¡” es la intensidad de lluvia efectiva, 
estimada como la intensidad de precipitación 
que no es interceptada por la vegetación, en 
el intervalo n de discretización temporal y 
9, el contenido volumétrico de humedad de 
saturación. El primer término de (5) representa 
la capacidad de infiltración y el segundo la 
disponibilidad de agua para infiltrarse. Se 
tiene prevista la desecación como condición 
de borde superior, aunque en la aplicación 
para este trabajo no fue considerada. Como 
condición de borde inferior, se consideran 
las celdas incluidas en la capa acuífera con 


humedad de saturación. El conjunto de celdas 
saturadas depende del nivel freático, y éste se 
actualiza dinámicamente en el modelo general. 

Para completar la solución del problema 
deben conocerse parámetros del suelo tales 
como las curvas de conductividad hidráulica 
versus humedad K(6), potencial mátrico versus 
humedad y(8) y difusividad versus humedad 
D(6). En el presente trabajo se han adoptado las 
relaciones de Brooks-Corey, representadas por 
las siguientes ecuaciones: 


A 
Ve ¡=> > 
9-0, (e) si => || = || (6) 


1 si[u|<|w 


donde 6, es el contenido de humedad co- 
rrespondiente a saturación residual; Á es un 
parámetro que depende de la distribución de 
los tamaños de poros; w, es la succión crítica 
(presión de entrada de aire), y S, es la saturación 
efectiva. Para la estimación de la conductividad 
hidráulica relativa, K_= K(0)/K, siendo K, la 
conductividad hidráulica saturada, se utiliza la 
ecuación propuesta por Corey (1977): 


uo el) (7) 


Las ecuaciones (6) y (7) se parametrizan 
para cada una de las celdas en que se divide el 
dominio de la ZNS. De esta manera se permite 
la representación de medios estratificados 
heterogéneos como los que se pretende 
estudiar. Los lazos de histéresis entre secado y 
mojadura no se tuvieron en cuenta. 


Funciones de pedotransferencia utilizadas 


Para estimar las propiedades hidráulicas de 
cada estrato se aplicaron funciones de pe- 
dotransferencia (FPT), constituidas por un 
conjunto de ecuaciones que relacionan pro- 
piedades hidráulicas de suelos con parámetros 
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físicos y granulométricos de las muestras de 
suelos. Existen numerosas FPT que actualmente 
se han incluido en códigos informáticos, como 
el SoilPar (Acutis éz Donatelli, 2003) y Rosetta 
(Schaap et al., 1999). En trabajos previos 
(Zimmermann é Basile, 2011) se utilizaron 
diferentes procedimientos para estimar 
propiedades hidráulicas de suelos limosos, 
partiendo de información granulométrica y 
de parámetros físico-químicos disponibles en 
las cartas de suelos de INTA. Se emplearon 
funciones de pedotransferencia (FPT) incluidas 
en los códigos SoilPar y Rosetta, junto con 
ecuaciones de regresión entre parámetros 
hidráulicos y parámetros físicos, dadas por 
un conjunto de autores. Se estimaron los 
parámetros hidráulicos medios ponderados 
para cuencas de la región y los pronósticos se 
compararon con valores obtenidos mediante 
calibración del modelo matemático hidrológico 
SHALL (Zimmermamn, 2003). 

A partir de los resultados obtenidos, pudo 
concluirse que, entre los códigos empleados, 
el que mejor se aproximó a los parámetros 
hidráulicos calibrados para la zona de estudio 
fue el programa Rosetta, suministrándole como 
información la composición textural, densidad 
aparente y dos puntos de la curva de retención 
(punto de marchitez permanente y capacidad 
de campo). A la vista de los resultados, se 
empleó entonces el código mencionado para 
la evaluación de los parámetros hidráulicos de 
cada estrato de los suelos analizados, lo que 
permite obtener las curvas de retención y de 
conductividad hidráulica mencionadas antes. 


Aplicación a suelos de la provincia de 
Santa Fe (pampa húmeda argentina) 


El área meridional de la provincia de Santa Fe 
(Argentina), motivo de interés en este estudio 
(figura 2), presenta un módulo pluviométrico 
del orden de los 1 000 mm y la temperatura 
media anual del orden de los 17 “C. Los mapas 
de suelos disponibles en la provincia de Santa 
Fe, elaborados por el INTA, clasifican, desde un 
punto de vista textural, los suelos de la región 


entre las fracciones franco limosa, franco arcillo 
limosa y arcillo limosa. De las cartas de suelos 
disponibles se extractó la información de 68 
perfiles de series tipo de suelos, representativos 
de los suelos presentes en el sur santafecino. 
En cada perfil de las series, en función de su 
profundidad y estratigrafía, se detallan entre 3 
y 8 horizontes, de los cuales se han publicado 
en dichas cartas la composición textural, 
junto con otros parámetros biofísicos, como 
contenido de materia orgánica, capacidad de 
intercambio catiónico, pH, etc. El conjunto de 
horizontes para todas las series estudiadas 
totalizan 434 muestras de suelo, considerando 
que, en promedio, cada serie de suelos está 
discretizada entre 6 y 7 horizontes. 


Estimación de parámetros por serie de suelos 


Con la información de las cartas pertenecientes 
a Cada perfil típico o serie de suelos, 
se confeccionaron planillas de cálculo, 
volcando características texturales y físico- 
químicas de cada estrato de las series. Con 
esta información, utilizando funciones de 
pedotransferencia intrínsecas del programa 
SoilPar, se obtuvieron densidades aparentes, 
humedades para capacidad de campo (CC), 
punto de marchitez permanente (PMP) y 
otras tensiones matriciales. Luego, se aplicó 
las FPTs incluidas en el programa Rosetta, 
obteniendo así las siguientes propiedades 
hidráulicas por horizonte o estrato de cada 
serie: conductividad hidráulica saturada, K; 
humedad de saturación, 6 ; humedad residual, 
9; conectividad de poros, A, y succión crítica, 
Y. 

Se ha considerado para la modelación de 
la ZNS un intervalo de discretización espacial 
Az de 0.05 m, lo cual da suficiente detalle para 
la estimación de perfiles de humedad en la 
región de análisis. Para la caracterización de 
las propiedades de cada celda se utilizó una 
rutina de interpolación lineal, que interpoló 
los parámetros ajustados para cada estrato. 
La zona que cubrió el dominio de modelación 
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Figura 2. Área de estudio: sector meridional de Santa Fe. 


alcanzó los 3.5 m de profundidad y, teniendo en 
cuenta el paso de discretización Az, dio lugar a 
70 celdas en sentido vertical, para cada una de 
las series de suelo estudiadas. 


Condiciones iniciales y de borde simuladas 
Las condiciones de borde e iniciales constituyen 


un factor importante en la simulación, porque 
los resultados obtenidos están influenciados 


por ellas. Como condición de borde superior 
se tomó la descrita anteriormente, esto es, 
condición de flujo preestablecido asociado con 
la capacidad de infiltración. Como condición 
de borde inferior se consideró el nivel freático 
a 3.5 m de profundidad. Esta profundidad, 
por un lado, no alterará en gran medida los 
fenómenos que se estudiarán en superficie y, 
por el otro, reflejará una situación frecuente en 
los sistemas hidrológicos estudiados. 
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La condición inicial se estableció teniendo 
en cuenta los siguientes criterios: a) en las pri- 
meras capas superficiales, hasta una profun- 
didad cercana a la extensión radicular de los 
cultivos presentes en los sistemas estudiados, 
la alteración del perfil de humedad se hace sig- 
nificativa y b) el perfil de humedad en profun- 
didad se ve escasamente alterado por las con- 
diciones ambientales externas. Como situación 
compatible con la de borde inferior fijada se 
estableció un perfil de humedad en equilibrio 
hasta el metro de profundidad (1 m aproxima- 
damente representa la profundidad radicular). 
Por encima del metro de profundidad, la va- 
riación del perfil de humedad se estableció en 
forma lineal (arbitraria y simplificada) entre el 
perfil de equilibrio y las diferentes condiciones 
de humedad prefijadas para el análisis: entre 
un 40 y 90% de la saturación efectiva (de 0.4 a 
0.9 de Se). La figura 3 muestra las condiciones 
adoptadas. 
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Figura 3. Perfiles de humedad para fijar las condiciones 


iniciales. 


Simulaciones realizadas 


Se realizaron simulaciones suponiendo cada 
perfil en condiciones de encharcamiento en 
superficie, logradas mediante una precipitación 
de elevada intensidad (70 mm/h) constante 
durante 12 horas de duración. Bajo estas 
condiciones se sometieron los 68 perfiles de 
suelos, analizando en cada caso la evolución de 
los perfiles de humedad en profundidad con el 
tiempo y, principalmente, la tasa de infiltración 
en el contorno superior. A manera de ejemplo, 
en las figuras 4 y 5 se presentan dos perfiles 
con características disímiles y comportamiento 
hidráulico consecuentemente diferenciado. 
Las características texturales de los mismos 
se presentan en el cuadro 1. Dichas figuras 
muestran los perfiles de humedad en diferentes 
tiempos de simulación, para las series “El 
Pedernal” y “Roldán” seleccionadas y para 
las condiciones extremas de humedad inicial. 
En conjunto se representan contenidos de 
humedad de saturación (8) y conductividades 
hidráulicas saturadas (K) en función de la 
profundidad (mediante interpolación lineal 
entre los valores asignados a cada estrato). 

La serie El Pedernal se encuentra en el 
sector sudoeste de la provincia, desarrollada 
a partir de sedimentos de textura franco 
limosa liviana, con importantes porcentajes 
de arena (cercanos al 18%). Esto le confiere 
una saturada 
relativamente elevada (cercana a 1 cm/h). 


conductividad hidráulica 


En la simulación, a las 48 horas, el frente de 
humedad alcanza profundidades superiores a 
los 2 m y su patrón es prácticamente uniforme. 
Los primeros estratos tienen mayor contenido 
de arcilla y retardan el flujo en profundidad, 
como puede observarse en la figura 4. La serie 
Roldán, de drenaje deficiente, se encuentra 
presente en el sector sureste de la provincia, en 
las cercanías de la ciudad de Rosario. Presenta 
un contenido de arcilla semejante (cercano al 
30%), distribuido en distintas profundidades. 
Las conductividades hidráulicas saturadas son 
del orden de 0.4 cm /h. En general, el alcance del 
frente de humedad a las 48 horas es del orden 
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Perfil tipo: el Pedernal 


Produndidad (m) 


0.2 0.4 0.6 0.8 1 1.2 
Ksat (cm /h) - Humedades (cm3 / cm3) 


Perfil tipo: el Pedernal 


Produndidad (m) 


Figura 4. Perfiles de humedad para distintos tiempos de simulación. Serie “El Pedernal”: a) Se inicial = 0.9; b) Se inicial = 0.4. 
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Perfil tipo: Roldán 


Produndidad (m) 


0.2 0.4 0.6 0.8 1 1.2 


Ksat (cm /h) - Humedades (cm3 /cm3) 


Perfil tipo: Roldán 


Produndidad (m) 


0.2 0.4 0.6 0.8 1 1.2 


Ksat (cm /h) - Humedades (cm3 /cm3) 


Figura 5. Perfiles de humedad para distintos tiempos de simulación. Serie “Roldán”: a) Se inicial = 0.9; b) Se inicial = 0.4. 
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Cuadro 1. Composición textural de los estratos para las series analizadas. 


Profundidades (m) Espesor (m) Arena (%) Arcilla (%) 
Serie El Pedernal 
<0.26 0.26 15.9 14.8 
0.26 0.40 0.14 12.2 17.3 
0.40 0.60 0.20 155) 7.6 
0.60 0.80 0.20 18.5 8.1 
0.80 115 0.35 22.6 92 
1.15 1.70 0.55 17.8 7.6 
% ponderado por espesor 17.6 919) 
Serie Roldán 
- 0.12 0.12 22 23.7 
0.12 0.25 0.13 2.4 26.3 
0.25 0.38 0.13 3.0 28.6 
0.38 0.80 0.42 2,2 50.7 
0.80 1.30 0.5 31) Sl 
1.30 1.75 0.45 3.3 19.8 
178 2.10 0.35 4.0 18.4 
% ponderado por espesor 3.1 30.1 


de 1 a 1.5 m, dependiendo de las condiciones 
iniciales de humedad. 


Abordaje simplificado para modelos de 
eventos 


El modelo SOLUM puede describir con detalle 
la distribución de contenidos de humedad 
del perfil de suelo y los flujos de intercambio 
entre la superficie y el nivel freático. En 
contraposición, requiere de un volumen 
apreciable de información y tiempo para 
ejecutar las simulaciones. Con frecuencia, en 
el empleo de modelos de eventos de lluvia 
no es necesario tanto detalle de resultados, 
sino que sólo con el conocimiento de los flujos 
de intercambio en la superficie del terreno 
(infiltración) es suficiente. La disyuntiva a 
resolver consiste en encontrar “parámetros 
equivalentes y únicos” de un suelo estratificado, 
que permitan estimar de forma razonable las 
tasas de infiltración de la superficie del terreno. 
El término “únicos” hace referencia a un medio 
equivalente homogéneo. De esta manera, 
podría aplicarse una metodología simple 


para evaluar la infiltración. En este trabajo, 
entonces, se ha propuesto el siguiente objetivo: 
sintetizar las heterogeneidades del suelo en 
valores representativos de las propiedades 
del conjunto. En particular, 
utilizando el modelo simplificado de frente 
húmedo (Green 8 Ampt), se trataría de obtener 
un valor de la conductividad hidráulica K, y 
del potencial de succión del frente húmedo y e 
ambos equivalentes a un suelo homogéneo, 


hidráulicas 


que produzcan el mismo volumen infiltrado 
que el suelo estratificado y heterogéneo. 


Método de Green y Ampt como algoritmo de 
pérdidas 


El método representa una simplificación del 
fenómeno real, que aproxima el avance de la 
humedad en el perfil de suelo mediante un fren- 
te húmedo, con condición de encharcamiento 
permanente en superficie. El suelo se considera 
homogéneo y profundo. Dicho frente es una 
) que 
divide el suelo con una humedad inicial, O, del 


AR 


frontera plana (flujo en forma de “pistón 


suelo saturado con un contenido de humedad 
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igual a la saturación, 6, y un potencial de 
succión y ¡en el frente de avance. 

Planteando para el movimiento vertical la 
ecuación de continuidad de masa y la ley de 
Darcy (generalizada para flujo no saturado) 
entre la superficie encharcada y el punto de 
avance del frente húmedo, se arriba a una 
formulación para la tasa de infiltración f y de 
la profundidad acumulada de agua infiltrada F 
en función del tiempo: 


TZ 
10 [to j (8) 
Flt 
F(t) =Kt+1p, AO afro 5 (9) 


donde K es la conductividad hidráulica vertical 
del suelo; y ñ el potencial de succión en el frente 
húmedo; A8, la diferencia entre contenidos de 
humedad de saturación e inicial A9 = 6. - 8, y 
t es el tiempo. La aplicación de este método 
requiere la estimación de la conductividad 
hidráulica K y el potencial de succión y, 


Obtención de los parámetros K. y Y,, 
equivalentes 


En principio, se propuso utilizar el modelo 
SOLUM para evaluar la tasa de infiltración 
en los estratos superficiales frente a una 
precipitación con intensidad y duración 
suficiente como para superar la capacidad de 
infiltración del conjunto. Si bien el modelo 
permite conocer la distribución de humedad en 
cada momento y a lo largo del perfil del suelo, la 
información de interés primordial que sintetiza 
el movimiento de humedad en el perfil es la 
curva de infiltración del primer estrato. Con 
la distribución temporal de la infiltración de la 
celda superficial f,(t) se ajustan los parámetros 
K, y , equivalentes para que la función [(t) 
evaluada con la expresión (8) y la f/(t) se 
aproximen de manera apropiada, utilizando 
un algoritmo de calibración propuesto en 
específico en este trabajo. 


Algoritmo de calibración 


Para la estimación de los parámetros K, y 
Y A equivalentes, fue necesario diseñar un 
algoritmo de ajuste, el cual fue incluido 
en el código computacional del programa 
SOLUM. Se propuso un algoritmo híbrido de 
calibración, que combina un “ajuste grueso” en 
el dominio restringido de los parámetros, junto 
con un algoritmo de “ajuste fino”, basado en el 
método del gradiente. 

La investigación preliminar se realiza 
mediante la evaluación de la función objetivo 
en una grilla trazada sobre la región de validez 
paramétrica subdividida en intervalos de cien 
partes. 

El algoritmo del gradiente para la búsqueda 
detallada de los óptimos parte de la definición 
de una Función Objetivo (FO) multidimensional 
(FO = y(x), donde x es el vector de parámetros). 
Básicamente, el algoritmo parte de un punto 
base posible ea 0 sobre el cual se evalúa 
la FO = y(x0), y el gradiente de la FO, Vy. Para 
construir el gradiente, se estiman las derivadas 
parciales mediante cocientes incrementales del 
tipo: 


Y _ il + Ax;,..X) 
0X; Ax; 


1 2 


TN 
Ax: 


1 


(10) 


para ¡=1,2,... n, donde Ax, es un valor pequeño 
respecto al valor de x,. El gradiente proporciona 
un vector que determina la dirección del 
óptimo de la FO. Se selecciona un nuevo punto 
en el espacio multidimensional cuya FO y(x') 
se acerca al óptimo buscado. En forma genérica 
puede expresarse que: 
x= x* + VyAh (11) 
El signo del segundo miembro se adopta 
positivo si la búsqueda es hacia un máximo 
y viceversa. El paso de búsqueda Ah, que 
representa la distancia al punto x', debe cumplir 
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con las restricciones impuestas alos parámetros. 
El proceso de búsqueda finaliza cuando en 
dos iteraciones sucesivas, la diferencia entre 
las FO calculadas es inferior a una tolerancia 
preestablecida y/o es superada una cantidad 
máxima deiteraciones fijada de antemano. En el 
caso particular de este estudio, los parámetros 
son dos (x,=Kx,=Y 1) y como función objetivo 
se consideró el error cuadrático medio (RSME) 
de las diferencias entre la tasa de infiltración 
en la celda superficial, f.(t), y la A(t) evaluada 
con la ecuación (8). Para la implementación del 
algoritmo se desarrolló con una subrutina en 
lenguaje FORTRAN, incorporada en el cuerpo 
principal del modelo SOLUM. 


Metodología de ponderación 


Paralelamente, se propuso una metodología de 
agregación para obtener parámetros equiva- 
lentes por serie de suelos, partiendo de los 
parámetros calculados para cada horizonte 
(Zimmermann et al., 2008). Dicha metodología 
es la que se pretende validar en este trabajo. 
Los criterios fueron los siguientes: a) para la 
conductividad hidráulica saturada vertical 
ponderada, K, se calculó la media armónica de 
las conductividades de cada estrato (ecuación 
(12) y b) para el potencial de succión del frente 
húmedo ponderado y, se adoptó el promedio 
ponderado por el espesor de cada estrato 
(ecuación (13)): 


L 
K,=3 (12) 
y 
¡=1 K; 
1 N 
Y, = [2% e; (13) 


donde K, y y, son los parámetros para los 
espesores €, N es el total de horizontes de la 
serie, y L es el espesor total, los cuales fueron 
obtenidos como se describió anteriormente. 
Para cada estrato, la conductividad hidráulica 
K,se consideró igual a la saturada y el potencial 


de succión para el frente húmedo, y, se estimó 
de la siguiente forma (Muñoz-Carpena éz 
Gowdish, 2005): 


Ñ 2 +3 Yo 


¿= (14) 
* 24+2 2 


donde la succión crítica de cada estrato y, la 
conectividad de poros A y la conductividad 
saturada K, se obtuvieron utilizando las FPT 
incluidas en el programa Rosetta, como se había 
mencionado. 


Comparativa entre ambas metodologías 


Los parámetros obtenidos mediante la meto- 
dología de ponderación propuesta (K, y w,) se 
compararon con los derivados de la calibración 
(K, y Ve equivalentes), previa modelación a 
través del programa SOLUM. 

Se presentó como alternativa en las com- 
paraciones integraciones de la ponderación, 
considerando distintas profundidades (L, en 
las ecuaciones (12) y (13)). En principio, se 
adoptó una profundidad de 2.5 m, conforme 
a la integración de la casi totalidad del perfil 
(considerado que la profundidad de la freática 
se propuso en 3.5 m en las simulaciones). 

Los valores de K, (en abscisas) y K, (en 
ordenadas) se grafican en la figura 6, para 
el conjunto de las series analizadas (68) y las 
distintas humedades iniciales consideradas 
e integradas hasta los 2.5 m de profundidad. 
Como se puede apreciar, los coeficientes de 
determinación obtenidos (e impresos en el 
ángulo superior izquierdo de cada figura) son 
relativamente bajos. 

En una segunda instancia fueron considera- 
dos los primeros estratos en las integraciones 
para la ponderación, aproximadamente los 
horizontes “A” y “B”, hasta una profundidad 
de 0.50 m. 

Los valores de K, (estimados) y K, (ajustados) 
se grafican en la figura 7, para las distintas 
series y humedades iniciales consideradas. 
Los coeficientes de determinación obtenidos 
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Figura 6. Comparaciones entre K, ponderados a una profundidad de 2.50 m (abscisas) y K, ajustados (ordenadas) para 
diferentes Se iniciales: a) 0.40, b) 0.50, c) 0.60, d) 0.70, e) 0.80 y £) 0.90, respectivamente. 
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Figura 7. Comparaciones entre Kp ponderados a una profundidad de 0.50 m (abscisas) y Ke ajustados (ordenadas) para 
diferentes Se iniciales: a) 0.40, b) 0.50, c) 0.60, d) 0.70, e) 0.80 y £) 0.90, respectivamente. 
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(e impresos en el ángulo superior izquierdo de 
cada figura) son considerablemente superiores 
a los anteriores. 

Esto último da la pauta que para los fines 
de caracterizar el proceso de infiltración en la 
superficie del terreno da mejor resultado, al 
considerar los primeros horizontes del suelo 
para parametrizar las propiedades hidráulicas 
del 
afirmarse que son los primeros estratos los 


conjunto. En otras palabras, podría 
que dominan el proceso de infiltración bajo la 
superficie. 

Para analizar con mayor detalle el efecto de 
las condiciones antecedentes de humedad, se 
obtuvo un valor promedio de conductividad 
hidráulica para las 68 series de suelos analizadas 
tanto para los valores ajustados como para los 
estimados mediante ponderación. Esto se hizo 
con el fin de contar con un valor representativo 
del conjunto de suelos analizados, aunque este 
promedio no tenga aplicación alguna salvo en 
el análisis propuesto en esta parte del trabajo. 

La figura 8 muestra la comparación de 
valores promedio ajustados y ponderados de 


conductividades hidráulicas para todas las 


series de suelo en función del contenido de 
humedad inicial. Allíse observa que los valores 
ajustados de K, promedio son asintóticos a los 
valores ponderados K, análogos para contenidos 
de humedad cercanos a la saturación. 

La ponderación inversa o armónica se utiliza 
para obtener conductividades equivalentes 
cuando el flujo circula en sentido transversal 
a un suelo saturado estratificado (Custodio dz 
Llamas, 1983). En el caso analizado, el medio 
no está saturado, por lo tanto los resultados 
se acercarán a la metodología de ponderación 
cuando más cerca se esté de las hipótesis 
de saturación. En la figura 8 se plasma un 
resultado lógico, que además permite analizar 
el grado de aceptación que se le puede brindar 
a la fórmula de ponderación. 

Por ejemplo, para las series de suelos 
analizadas, serían aceptables las estimaciones 
de conductividades hidráulicas saturadas 
mediante ponderación hasta un rango del 
orden del 50% de su saturación efectiva con un 
rango de errores acotados al 10%. 

Si se piensa, además, que en la región 
analizada, los tenores de humedad para las 
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Figura 8. Comparación de valores promedio de Kp y Ke para todas las series de suelos en función de la humedad inicial. 
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profundidades hasta el horizonte “B” en la 
generalidad de los casos se encuentran por 
encima del 50% de la saturación efectiva 
debido al carácter “húmedo” del clima, 
puede aseverarse que la metodología de 
ponderación es adecuada para representar las 
conductividades hidráulicas del conjunto de 


estratos. 
Conclusiones 


Se ha validado una metodología de agregación 
para obtener conductividades hidráulicas equi- 
valentes de un medio homogéneo en suelos 
heterogéneos no saturados. El procedimiento 
consistió en dos etapas. En una primera etapa 
se propuso una modelación numérica del flujo 
en zona no saturada, específicamente para este 
trabajo, en el dominio heterogéneo discretizado 
con alto grado de detalle y parametrizado 
mediante funciones de pedotransferencia 
adecuadas. En una segunda etapa se calibraron 
los parámetros equivalentes del modelo de 
frente húmedo con los resultados alcanzados 
en la primera etapa, comparándolos con los 
propuestos en la metodología de agregación. 
El procedimiento se aplicó a 68 perfiles típicos 
de suelos que caracterizan el sector meridional 
de la provincia de Santa Fe (Argentina). Los 
resultados alcanzados validan la metodología 
de agregación propuesta para un numeroso 
muestreo de composiciones texturales, pro- 
piedades hidráulicas y contenidos de humedad 
antecedente. Con el fin de caracterizar el 
proceso de infiltración en la superficie del 
terreno, da mejor resultado considerar en 
las ponderaciones los primeros horizontes 
del suelo para parametrizar las propiedades 
hidráulicas del conjunto. Podría afirmarse que 
son los primeros estratos los que dominan el 
proceso de infiltración bajo la superficie. 

Se observa que los valores calibrados de con- 
ductividades hidráulicas son asintóticos a los 
valores ponderados análogos para humedades 
cercanas a la saturación. Considerando que 
la ponderación armónica propuesta se utiliza 
para obtener conductividades equivalentes 


cuando el flujo circula en sentido transversal a 
un suelo saturado estratificado, los resultados 
serán mejores cuando más cerca se esté de las 
hipótesis de saturación. Para las series de suelos 
analizadas serían aceptables las estimaciones 


de conductividades hidráulicas saturadas 


mediante ponderación hasta saturaciones 


efectivas del orden del 50%, con un rango de 
errores acotados al 10%. 
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